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A  mini conference  on  the  topic  of  Order  Restricted  Inference 
was  held  on  the  University  of  Iowa  campus  during  the  days  of 
September  11,  12,  and  13,  1983.  The  conference  was  supported  by 
a  grant  from  the  Office  of  Naval  Research  and  the  University  of 
Iowa.  The  conference  consisted  of  twenty  presentations  on  a  wide 
variety  of  topics  dealing  with  Order  Restricted  Inference.  There 
were  approximatel y  35  conference  participants  coming  from  a  wide 
variety  of  backgrounds  and  locations. 


7=^  The  majority  of  papers  presented  at  the  conference  are 
currently  being  organized  to  appear  in  a  Conference  Proceedings 
Volume  to  be  published  by  Springer — Verlag  during  the  summer  of 
1986.  ^ 
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CONFERENCE  SUMMARY 


The  Office  of  Naval  Research  issued  Grant  No.  N00014-G-0161 
for  the  purpose  of  holding  a  miniconference  on  the  topic  of  order 
restricted  inference.  This  conference  was  held  in  the  Iowa 
Memorial  Union  on  the  University  of  Iowa  campus  on  September  11, 

12,  and  13,  1985.  There  were  approximately  thirty-five  participants 
at  the  conference  and  the  program  consisted  of  twenty  presentations. 
Enclosed  is  a  list  of  participants  for  the  conference,  a  copy  of 
the  program  of  the  conference,  and  a  partial  list  of  abstracts  for 
the  presentations  given.  Funds  from  the  Office  of  Naval  Research 
grant  were  used  primarily  to  support  travel  and  subsistence  of 
twenty-one  of  the  participants. 

We  feel  the  conference  was  an  unqualified  success  and  we  have 
received  numerous  favorable  comments  from  the  participants.  Many 
of  the  participants  indicated  that  if  a  similar  such  conference 
is  to  be  held  in  the  future,  they  would  wish  to  attend.  (A 
similar  conference  on  the  same  topic  was  held  at  the  University  of 
Iowa  in  the  spring  of  1981  with  support  from  ONR. )  The  tone  of 
the  conference  was  enthusiastic,  and  the  interaction  between 
participants  was  outstanding. 

It  has  been  decided  to  publish  a  Conference  Proceedings  which 
will  include  many  of  the  best  papers  of  the  conference.  With  this 
in  mind,  we  contacted  several  publishing  concerns  incuding  Marcel 
Dekker,  Springer-Verlag,  Academic  Press,  and  the  IMS  Lecture  Note 
Series  to  see  if  they  would  have  an  interest  in  such  a  publication. 
The  response  of  the  publishing  concerns  was  positive,  and  we  have 
decided  to  accept  an  offer  by  Springer — Verlag  to  produce  such  a 
Proceedings.  We  anticipate  that  this  volume  will  be  available 
during  the  late  summer  of  1986. 
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Testing  Whether  a  Set  of  Normal  Means  Satisy  a  Partial  Order 

Roger  L.  Berger 

North  Carolina  State  University 


Let  denote  l  normal  means.  Robertson  and  Wegman 

(1S78,  Annals  of  Statistics)  consider  testing  the  null  hypothesis 
Hj:  <vi  >  iE  isotone  with  respect  to  a  specified  partial 

order  versus  the  alternative  H  2 :  not  H^.  In  this  paper  the 
likelihood  ratio  test  (LPT)  for  the  problem  of  testing  the  null 
hypothesis  H2  versus  the  alternative  Hj  is  discussed.  If  one 
wishes  to  particularly  guard  against  deciding  the  means  are 
ordered,  when  they  are  not,  our  formulation  may  be  preferable  to 
Robertson  and  Wegman’ s.  The  LRT  we  obtain  is  quite  different  from 
Robertson  and  Wegman ’s.  Our  test  statistic  is  a  minimum  of  t- 
steti sties  and  the  critical  values  for  our  test  are  central  t- 
distribution  percentiles.  Also  the  least  favorable  configuration 
in  our  problem  is  not  the  configuration  in  which  all  the  means 
are  eoual . 


Testing  Homogeneity  of  Means  Against  The  Alternative  of  Simple 

Order  From  A  Bayesian  Viewpoint 

Bahadur  Singh 
Memphis  State  University 


nl  for  i  =  1,2,. 


Assume  the  CY..  ;j=l,  ... 
l] 

independent  random  samples  from  the  normally  distributed 


k  are  k 


population  with  mean  y.  and  common  variance  o  .  It  is  desired 

to  test  the  null  hyposthesis  Ho:  p  j  =  p2  =  •••  =  P^  against 

the  alternative  hypothesis  of  simple  order  H^:  pi<p2  <•.. 

<  with  at  least  one  strict  inequality.  Letting  8^=  P^+j  ~  y^ 

(i  ~  1,  2 ,  ...,  d-1 ) ,  and  assuming  independent,  noninf ormati ve 

prior  distributions  on  the  unknown  parameters,  the  postierior 

distribution  of  0  =  (  ex,  8^.  ^  )  is  obtained  for  the  two 

2 

cases  of  known  and  unknown  o  .  In  the  notation  described  above 

testing  H  versus  Hi  -H  is  equivalent  to  testing  the  null 
o  o 


hypostesis  K  :  e  i  =  6  =  .  .  .  =  e 


k-1 


=  0  versus  the  alternati 


hyposthesis  K 


0  ^  >  0  for  some  <i  =  1,2,...,  k-1).  An 


approximate  and  asymptotic  chisquare  test  statistic  is  derived 

2 

for  testing  K  versus  K  i  when  is  o  known,  while  an 
o  * 

approximate  and  asymptotic  F  test  statistic  is  derived  for 


testing  k  versus  K  ■  when  a  is  unknown.  The  Bayesian  test 
o  1 

procedures  are  developed  using  the  concept  of  the  highest 


posterior  desity  regions. 


“2 

Robustness  of  % 
Respect 


~z 

-  Like  and  E  -  Like  Statistics  with 
to  Departures  from  Normality 

Michael  J.  Schell  and  Bahadur  Singh 
St  Jude  Children's  Research  Hospital 

Let  i  v..:  i=l . k;  j=l,...,n>  be  k  independent  random 

13  .  2 
samples  of  size  n  with  means  y^  and  common  variance  o  ■  It  is 

desired  to  have  tests  for  Hg  :  y  i  -  p  2  =‘  ■  ■  ~  *  versus  the 

alternative  of  simple  order  Hj  s  y^<...  <11^ »  with  at  least  one 

strict  inequality  and  Hj  versus  H2 ‘  y^>y^_^  for  some  i.  When 

the  k  populations  are  normal,  the  likelihood  ratio  statistics 

2  2 
are  T0 j  and  Tj2  when  0  is  known  and  Sg *  and  S  12  when  o  is 

unknown.  Tgj  and  Sg j  are  commonly  known  as  chi  — bar-square  and 

E-bat — square,  respectively.  Using  simulations  methods,  this 

paper  assesses  the  robustness  of  the  significance  levels  of  these 

tests  wher.  the  k  papulations  are  non-normal  .  The  distributions 

used  are  the  uniform,  the  double  exponential,  the  t  with  3 

degrees  of  freedom  and  mixtures  of  normals.  In  addition, 

estimates  for  the  power  functions  of  the  selected  non-normal 

distributions  are  compared  to  the  corresponding  power  function 

based  on  an  underlying  normal  distribution. 


Dose-Response  Analysis  Under 
Unimodalit>  of  Response-to-Dose 

I 

Richard  L.  Schmoyer 
Oak  Ridge  National  Laboratory 

A  dose-response  curve  it  is  sigmoidal  if  it  is  nondecreasing  and 
there  is  a  point  M  to  the  left  of  which  it  is  convex  and  to  the 
right  of  which  ir  is  concave.  If  tt  is  sigmoidal,  then  the 
response-to-dose  ratio  ir(x>/  x  is  uni  modal  and  maximized  at  a 
point  H  JiM.  If  *  (0)  is  known,  if  H  is  known,  and  if  xQ  x  «£  H, 
then  an  upper  confidence  bound  (UCB)  for  tt  (  Xg)  can  be  constructed 
from  a  UCB  far  n  (x).  Unfortunately,  H  is  seldom  known  in 
practice.  In  the  setting  of  quantal  bioassay,  two  lower 
confidence  bounds  for  H  and  correspondi ng  UCB’s  for  it  (Xq  )  when  H 
is  unknown  are  obtained  and  contrasted.  One  is  based  on  a 
likelihood-ratio  statistic;  the  other  on  a  weighted  sum  of 
constraint  contrasts.  A  maximum-violator  algorithm  with 
guaranteed  convergence  is  given  for  computing  the  maximum 
likelihood  estimates  of  the  response  probabilities  subject  to 
unimodality  of  it  (x)/  x. 

Keywords:  accelerated  testing;  linear  i nterpol at i on ;  quantal 

bioassay. 


Conditional  Tests  with  an  Order 
Restriction  as  a  Null  Hypothesis 

Peter  C.  Moll  an  and  Richard  L.  Dykstra 
Michigan  Technological  University  &  University  of  Iowa 

For  the  isotonic  normal  means  problem,  Bartholomew  (1961) 

disc  tssed  a  conditional  1 i kel ihood-rati o  test  of  H  :  the  means 

o 

are  homogeneous,,  vs.  H  :  the  means  satisy  the  linear  order.  He 

1 

concluded  that  the  conditional  test  was  substantialy  less 
powerful  than  the  chi-bar — squared  test.  However,  for  testing  H^ 
vs.  H2:  all  alternatives,  the  correspond! ng  conditional  test  can 
be  more  powerful  than  the  chi -bar — squared  test.  Moreover,  the 
conditional  test  can  be  modified  so  as  to  be  asymptotically 
similar. 

1 

These  conditional  tests  are  of  particular  interest  in  general 
tests  of  simultaneous  inequality  constraints  on  parameters  of 
asymptotically  normal  distributions,  for  which  the  coefficients 
corresponding  to  the  p(l,k)’s  are  difficult  to  obtain.  In  this 
general  context,  the  likelihood  ratio  statistic  is  asymptotical  1 y 
chi -bar — squared  whenever  the  true  parameter  vector  lies  in  H-^;  we 
outline  a  new  proof  based  on  Si  Ivey’s  theorem  that  a  constrained 
estimate  and  its  correspondi ng  vector  of  Lagrange  multipliers  are 
asympototical ly  normal  and  independent. 


Statistical  Inference  Under 
Complete  Monotonicity  Restrictions 


Douglas  R.  Miller  and  Ariel  a  Sofer 
The  George  Washington  Univ.  &  George  Mason  Univ. 

A  completely  monotone  -function  is  a  -function  which  has 
nonnegative  even-numbered  derivatives,  and  nonpositive  odd- 
numbered  derivatives.  In  this  paper  we  address  the  problem  o-f 
regression  under  the  restrictions  of  complete  monotonicity. 

These  problems  may  arise  in  various  applications  and  one  such 
application  in  software  reliability  is  described. 

« 

The  paper  discusses  some  of  the  numerical  difficulties  associated 
with  these  problems.  In  constrast  to  isotone  regression,  where 
the  pool  adjacent  violators  algorithm  solves  the  problem  in  a 
finite  number  of  simple  steps,  there  is  no  simple  finite 
algorithm  for  solving  weighted  least  squares  when  the  additional 
constraints  of  isotone  (or  antitone)  higher  derivatives  (up  to 
some  maximal  order)  are  imposed.  We  show,  that,  as  the  problem 
size  and  the  maximal  order  of  derivatives  considered  grow,  the 
problem  becomes  increasingly  ill-conditioned.  The  various 
difficulties  which  may  arise  are  presented,  and  the  performance 
of  various  numerical  algorithms  discussed. 


Power  Series  Approximations  to  the  Null 
Distributions  in  Order  Inference: 

The  case  of  Unequal  Weights 

Bahadur  Singh  and  F.  T.  Wright 
Memphis  State  University  &  University  of  Missouri-Rol la 

Bartholomew’s  statistics  for  testing  homogeneity  of  normal  means 
with  ordered  alternatives  have  null  distributions  which  are 
mixtures  of  chi -squared  or  beta  distributions  depending  on 
whether  the  variances  are  known  or  not.  If  the  ratio  of  sample 
sizes  to  the  variances  are  not  equal,  the  mixing  coefficients 
can  be  difficult  if  not  impossible  to  compute.  Approximations  to 
the  significance  levels  of  these  tests  have  been  developed  for 
the  total  order  and  simple  tree  restrictions.  However,  for  a 
moderate  or  large  number  of  means,  these  approximations  can  be 

t 

tedious  to  implement.  Two  and  four — moment  expansions  in  terms 
of  Laguerre  and  Jacobi  polynomials  are  developed  to  facilitate 
the  use  of  these  approximations.  Approximate  significance  levels 
are  also  developed  for  the  testing  situation  in  which  the  order 
restriction  is  the  null  hypothesis.  Numerical  studies  show  that 
in  each  of  the  cases  the  two-moment  approximation  is  quite 
satisfactory. 


Bayes  and  Maximum  Likelihood  Order— Restri cted  Inference  for 
Models  for  Ordinal  Categorical  Data 


Alan  Agresti 

Department  of  Statistics 
University  of  Florida 

A  class  of  association  models  for  contingency  tables  has 
parameters  that  can  be  interpreted  as  category  scores.  For 
classifications  having  ordered  categories,  it  is  often  reasonable 
to  assume  that  the  score  parameters  have  a  corresponding 
ordering.  This  article  proposes  order — restricted  estimates  of 
score  parameters  in  these  models.  For  these  estimates,  the  local 
log  odds  ratios  have  uniform  sign.  Two  solutions  are  given.  For 
the  maximum  likelihood  solution,  goodness-of-f it  statistics  are 
•"®lated  to  statistics  for  collapsed  tables  and  to  statistics  for 
testing  equality  of  score  parameters.  For  the  Bayes  solution, 
prior  distributions  induce  the  order  restriction,  and  prior 
beliefs  reflecting  strong  association  have  the  effect  of  moving 
the  estimates  away  from  the  boundary  of  the  restricted  parameter 


space. 


Imputing  Missing  Categorical  Data 
under  the  Assumption  of  Positive  Association 

Jon  H.  Lemke  and  Gregg  A.  Drube 
University  of  Iowa 

Imputation  classes  are  subpopulations  of  a  surveyed  population 
for  which  one  is  willing  to  make  the  assumptions  that  item 
nonresponse  occurs  at  random,  even  though  the  data  cannot  be 
assumed  to  be  missing  at  random  for  the  entire  population.  For 
the  cross— classification  of  ordinal  categorial  variables,  we 
demonstrate  how  the  iterati ve-incremental  algorithm  of  Lemke  and 
Dykstra  <1984)  can  be  used  as  the  M-step  of  the  EM-algorithm  to 
obtain  maximum  likelihood  estimates  of  the  cell  probabilities 
within  imputation  classes  under  multiple  closed  convex  cone 
restrictions.  These  order  restrictions  can  be  either  on  the 
logarithms  of  the  cell  probabilities  or  simply  on  the  cell 
probabilities.  The  specific  application  for  this  presentation  is 
the  restriction  that  the  local  odds  ratios  are  all  at  least  one 
for  the  cross-classification  of  two  ordinal  response  variables, 
that  is,  we  are  assuming  that  within  each  imputation  class  there 
is  a  positive  association  between  the  two  response  variables. 


’Is- 


Shrinkage  Estimates  in  Isotonic  Regression 
A.  G.  Warrack 


When  estimating  a  set  of  means  from  k  independent  normal 
populations,  the  techniques  of  isotonic  regression  may  be  used  if 
it  is  known  that  the  means  satisfy  a  partial  order.  If  the 
partial  order  happens  to  be  a  simple  order,  the  isotonic  estimate 
of  the  means  has  the  desirable  property  that  its  expected  mean 
square  error  is  less  than  the  unrestricted  maximum  likelihood 
estimate,  and  this  property  holds  not  only  overall,  i.e.  when 
summed  over  all  k  populations,  but  also  pointwise,  for  each  of 
the  k  estimates.  ' 

However,  the  isotonic  estimates,  particulari ly  those  of  the 
largest  and  smallest  population  means,  are  biased.  In  an  effort 
to  reduce  this  bias,  shrinkage  techniques,  as  developed  by  Stein 
and  James,  are  investigated,  and  various  "shrunken”  isotonic 
estimates  are  compared  with  the  regular  isotonic  estimates. 

While  these  estimates  do  not  remove  the  bias  from  the  regular 
isotomic  estimates,  Monte  Carlo  results  are  presented  to  indicate 
that  the  overall  mean  square  error  of  the  isotomic  estimates  can 
be  reduced  by  shrinkage  techniques,  in  some  cases  considerably 

so. 
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An  asymptotical  1 y  di stri but i on— f ree  test 
for  ordered  alternatives  in  two-way  layouts 

Z.  Govi ndaraj ul u  and  S.H.  Mansouri -Ghiassi 
University  of  Kentucky 

An  asymptotically  distribution-free  test  is  proposed  for 
unequally  spaced  ordered  alternatives  in  two-way  layouts.  The 
test  statistic  is  a  linear  function  of  the  ranks  of  residuals 
when  the  nuisance  parameters  are  estimated.  We  show  that  the 
limiting  distribution  of  the  test  statistic,  when  properly 
standardized,  is  normal .  The  asymptotic  relative  efficiency 
comparisons  (in  Pitman  sense)  with  respect  to  the  likelihood 
derivative  test  and  nonparametric  tests  for  randomized  complete 
blocks  show  that  our  procedure  is  generally  more  powerful. 


An  Approach  to  Fitting  Convex  Interpolating 
and  Smoothing  Cubic  Splines 

Pat  Boyle  and  Richard  Dykstra 
University  of  Iowa 

Often  it  is  desirable  to  find  solutions  to  infinite  dimensional 

optimization  problems  which  have  an  infinite  number  of 

\ 

constraints.  An  algorithm  developed  by  Dykstra  (19B3)  is 
described  and  shown  to  be  useful  in  determining  convex 
interpolating  and  smoothing  functions  as  solutions  to  such 
problems.  Some  numerical  results  are  presented  based  upon  APL 
routines  that  were  written  to  obtain  the  convex  smoothing 
functions.  Finally,  some  concluding  observations  and  conjectures 


are  offered 


Applications  of  Isotonic  Regression 
*  £o  Multivariate  Density  Estimation 

Thomas  W.  Sager 

Center  for  Statistical  Sciences 
The  University  of  Texas  at  Austin 

Unrestricted  nonparametr i c  multivariate  density  estimation 
suffers  from  difficult  convergence  and  computational  problems. 
One  way  to  overcome  these  problems  is  to  exploit  presumed  or 
^estimated  structure  in  the  density.  The  isopleth  density 
^estimator  presumes  or  estimates  the  structure  of  the  contours 
®f  the  . density  to  effectively  reduce  the  dimensionality.  The 
.  -estimator  -incorporates  an  order— preserving  algorithm  to  insure 

^4that  tiigher  isopleths  ■have  higher  density  .estimates  than  lower 

/'  ;  .  "  .  '  -  .  ■  •  ;  *  if 

isopleths.  Convergence  properties  and  stimulations  are 
.presented.  The  importance -of  ,-«dge  effects  is  also -noted. 


Isotonic  M-Estimation 
Sue  Leurgans 

The  Ohio  State  University 

The  isotonic  regression  of  Y  on  x  is  the  isotonic  function  of  ;< 
which  minimizes  the  error  sum  of  squares.  For  linear  models, 
studies  of  robustness  properties  have  shown  that  the  use  of  other 
criterion  functions  is  valuable  in  practice.  One  large  class  of 
other  alternative  criterion  functions  gives  M-esti mators.  Magel 
and  Wright  <1984)  demonstrate  that  analogous  benefits  can  be 
obtained  for  isotonic  regression.  In  this  paper,  several 
recursi ve-partitioning  algorithms  which  are  valid  for  the  regular 
least  squares  problems  do  not  apply  to  all  M-estimation  problems, 
even  those  based  on  estimators  with  monotone  influence  functions. 
Through  examples,  isotonic  M-estimators  are  shown  to  implicitly 
balance  the  order  information  supplied  by  a  data  point  against 
the  size  of  the  residual  generated  by  the  point. 
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A  miniconference  on  the  topic  of  Order  Restricted  Inference  was 
held  on  the  University  of  Iowa  campus  during  the  days  of 
September  11,  12,  and  13,  1985.  The  conference  was  supported  by 
a  grant  from  the  Office  of  Naval  Research  and  the  University  of 
Iowa.  The  conference  consisted  of  twenty  presentations  on  a 
wide  variety  of  topics  dealing  with  Order  Restricted  Inference. 
There  were  approximately  35  conference  participants  coming  from 
a  wide  variety  of  ba 
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The  majority  of  papers  presentee  at  the  conference  a 
currently  being  organized  tc  appear  ir.  a  Conference 
Proceedings  Volume  to  be  pubiisnec  by  Sprmger-Veria 
during  the  summer  of  198c. 
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